The Construction of Sudden Cosmological Singularities 
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Solutions of the Friedmann-Lemaitre cosmological equations of general relativity have been found 
with finite-time singularities that are everywhere regular, have regular Hubble expansion rate, and 
obey the strong-energy conditions but possess pressure and acceleration singularities at finite time 
that are not associated with geodesic incompleteness. We show how these solutions with sudden 
singularities can be constructed using fractional series methods and find the limiting form of the 
equation of state on approach to the singularity. 
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, In standard Friedmann-Lemaitre cosmological models a fluid is placed in a homogeneous and isotropic spatial geom- 

■ etry whose dynamics is then determined by two independent Einstein equations for three unknown time-dependent 
H functions, the Friedmann metric scale factor, a{t), the fluid density, p{t), and fluid pressure, P{t), respectively (with 

units chosen with c = SttG = 1): 

a2 3 ^ ' 

O; _2«-^ = P. (2) 

^ If an equation of state P = f(p) is chosen, the system closes and the two remaining unknown functions are determined. 
Physically reasonable equations of state, like those for perfect fluids, produce well behaved expansion factors with 
behaviours that offer simple Newtonian interpretations. However, the challenge of providing a compelling explanation 
for the observed acceleration of the universe in terms of a 'dark energy' fluid has led to an exploration of other 
K*" , less familiar equations of state pj], motivated by the form of bulk viscous stresses 0, situations were dP/dp is 
not everywhere continuous, or scenarios in which there is no equation of state at all. These introduce quite different 
possibilities into the solution space of Friedmann universes and can produce unexpected types of finite-time singularity 
which are far softer than the curvature singularities studied previously (sj in connection with geodesic incompleteness 

Barrow [a, |7| identified a whole new class of pressure-driven singularities that keep the scale factor, a, expansion 

■ rate, a/a, and the density, p, finite while the pressure, P, blows up at a finite time despite the energy conditions 

\ p > and p + 3P > holding. Their status as stable solutions of the classical Einstein equations in the presence of 
small scalar, vector and tensor perturbations has been studied in a gauge covariant formalism |8|], and they have also 
been found to be stable against quantum particle production processes |9|] . They have been investigated in a number 
of different cosmological scenarios and their behaviour has been classified in the light of other types of finite-time 
5^ ' singularity that can arise in isotropic and anisotropic cosmologies ,lft-il4i]. The formal definitions of these singularities 
has been discussed by Lake [l^ in the case where infinities occur in second derivatives of a{t), but similar examples 
exist where the sudden singularity occurs in higher derivatives of a{t) and no energy conditions are threatened [3,ll3|- 
Here, we shall provide an alternative systematic way of constructing them. 

A sudden singularity will be said to arise everywhere at comoving proper time ts in a Friedmann universe expanding 
with scale factor a{t) if 

lim a{t) = =/= 0, lim d{t) = ds < oo, lim d(t) = oo, 

for some ts > 0. From the last requirement we have that, for all constant Af > 0, there exists a constant e > such 
that for |t — ts| < e and t < tg, we have 

d{t) > M. (3) 

By integrating the differential inequality ([3]), we have 



0<^ Mdt < ddt = ds-d{t), 
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and with a further integration we see that for every positive number M we can find a sufficiently small left neigh- 
bourhood of tg such that for every t in that neighborhood the scale factor satisfies 

a{t) > as" hslts -t) + Yi'ts- tf- 

For this condition to be satisfied in the neighborhood of ts, we must have 

a(t) = as- asits - t) + Cits - t)" + ■■■ (4) 

where C is a constant that depends on M, and necessarily we have 1 < rt < 2. 

If 71 is a rational, the expansion factor can be written as a Puiseux series of the form 



i=0 



where s > is a natural number. 

We can study these kinds of singularities by the dominant balance method p^ . wherein a dynamical system 
X = f(x) is split into dominant and subdominant terms f(x) = f(°)(x) + f '^^"'^^ (x) and then we search for solutions of 
the dominant dynamical system x = f''°)(x) in the form of a Puiseux or ^I* series. In our case, setting x = a, y = a 
we get from ^ 

k + y^ + x^P 

x = y, y^ — , (5) 

and we require a dominant solution x = and y — hs both finite, with pressure being infinite. Clearly, these two 
conditions cannot be simultaneously satisfied. 

The root of the problem is that the dominant balance method looks for a singularity of the variable 

^ ^ [ai{ts - tY\a2{ts - tY\ - ■ 

(i.e. at least one of the pi < 0) but we want x to be finite and x (and so d) to have a singularity since the pressure 
term appears in the higher derivative (o) equation Despite this structure, the general form of the pressure which 
leads to this singularity can be found by looking for solutions of ([S]) of the form 

oo oo 

s(t) = ^CH (^,-^r/^ y(i) = ^c2, (^«-^r/^ (6) 

with s > 1. We want 

cio = a 7^ 0, C20 = f3 ^ 0, 
to hold. In this case, the pressure takes the form 

oo 

p = {ts-tYJ2p^{ts-ty/% (7) 

i=0 

and for P to blow up at tg, we must have h < 0. The recursion relations for the Cy are 

s + i 

Cis = 0, 0<1<S, C2i ^ Ci(s+i), 



and 



2(s + l)a 



"Cl(s+ 



4(s + 2)a ,N-(s-2)/s 

l)(ts~t) Ci(s+2)(ls - 



^2(£_4_1) + - «(2s + l)ci(2.+i)](is - t)'/^ + 

[2(5 + 2)2/3 (s + i)(s + 3) 2 4(s + l)(s + 2)a 

+ ^2 ^2 ^l{s+l) ^2 Ci(2s+2) 

+ ■ • • = a^poits - + a^Piits - + c^P2{ts - tf+^'' + ■■■ 



[ts " tf'' 
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By balancing the terms in this equation we can determine the series expansion of the pressure that wih lead to a 
specific sudden singularity. For example, if n = 3/2 in Q 



we get for the pressure: 

P 
with 

and the density has the form 



a = as, Oi = (os - 1) ( ^ ) ( ) , ^ = -ai, 



i \ t 



We note here that the asymptotic form of the equation of state has the simple form 

lnTi(P(p-p.)) = ^, 

and so the fluid behaves like a Chaplygin gas on approach to the sudden singularity. Hence, we see that these 
formal methods for the analysis of finite-time singularities in dynamical systems allow us to construct cosmological 
singularities of the sudden type. 

In conclusion, we have shown how to construct evolutionary behaviours meeting the conditions for a sudden cos- 
mological singularity using analytic techniques with fractional power series employed to study finite-time singularities 
in differential equations. 
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